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ABSTRACT 



Context. The polarization sensitivity of the upcoming millimetric observatories will open new possibilities for studying the properties of 
galaxy clusters and for using them as powerful cosmological probes. For this reason it is necessary to investigate in detail the characteristics 
of the polarization signals produced by their highly ionized intra-cluster medium (ICM). This work is focussed on the polarization effect 
induced by the ICM bulk motions, the so-called kpSZ signal, which has an amplitude proportional to the optical depth and to the square of the 
, tangential velocity. 

f"*) ■ Aims. In particular we study how this polarization signal is affected by the internal dynamics of galaxy clusters and what is its dependence on 
the physical modelling adopted to describe the baryonic component. 

Methods. This is done by producing realistic kpSZ maps starting from the outputs of two different sets of high-resolution hydrodynamical 
N-body simulations. The first set (17 objects) follows only non-radiative hydrodynamics, while for each of 9 objects of the second set we 
(~~^ implement four different kinds of physical processes. 

• • ■ Results. Our results shows that the kpSZ signal turns out to be a very sensitive probe of the dynamical status of galaxy clusters. We find that 
e ^ ' major merger events can amplify the signal up to one order of magnitude with respect to relaxed clusters, reaching amplitude up to about 100 
k>( , nK. This result implies that the internal ICM dynamics must be taken into account when evaluating this signal because simplicistic models, 
j_J ■ based on spherical rigid bodies, may provide wrong estimates. In particular, the selection of enough relaxed clusters seems to be fundamental 
Cd to obtain a robust measurement of the intrinsic quadrupole of the cosmic microwave background through polarization. Finally we find that the 
dependence on the physical modelling of the baryonic component is relevant only in the very inner regions of clusters. 



> 

o 



O 



Key words, polarization - cosmic microwave background - galaxies: clusters - cosmology: theory - methods: numerical 



1. Introduction 

The last generation of X-ray satellites like Chandra and XMM- 
Newton have shown a complex scenario for the intracluster 
medium (ICM). In particular the presence of col d fronts has 



been detected in many galaxy clusters (see, e.g.. lAllen et al. 



200ll) . The existence of this phenomenon suggests that also 
relaxed clusters display a picture in contrast with the previ- 
ously hypothesized hydrostatic state for the ICM, showing mi- 
nor merger events and evident gas motions. 

In this paper we focus our attention on a different signal 
which could give strong constraints on the merging activity 
of galaxy clusters: the kine matic Sunyaev-Zel'dovich (SZ) po- 
larization (hereafter kpSZ) dSunvaev & Zerdovichlll980bl) . It 
takes origin by single scattering of moving electrons, and its 



dependence on the square of their tangential velocity makes 
it a suitable probe for investigating the dynamic properties of 
clusters during the non-linear process of their formation. 

As we will discuss in more detail in Sect. [2] galaxy clus- 
ters produce other linear polarization contributions to the cos- 
mic microwave bac kground (CMB) radiation thr ough different 



mic microwave background (ClVLhs ) radiation througn dirterent 
processes (see, e.g.. ISazonov & Sunyaev . 1999). Basically all 
components of the quadrupole moment present in the incom- 
ing radiation induce a linear polarization: in fact there is the 
polarization induced by the intrinsic quadrupole of the CMB 
radiation (Pq), by the quadrupole component of kinetic and 
thermal SZ effects {Ptp^-sz an d P tp n -sz, respectively) and by the 
gravitational lensing effect (P/) of clusters them-self. These po- 
larization contributions, with the exception of Pk P 2 sz> which 
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is one order of magnitude smaller, h ave amplitudes similar to 
the kpSZ effect dLavaux et all 120041) . Thus, none of them can 
be neglected. In any case, a component separation is possible, 
thanks to their different frequency dependences. It is important 
to notice that each of these signals carries different and impor- 
tant information on galaxy clusters and cosmology. For exam- 
ple, Pq could be used to estimate the intrinsic quadrupole at 
different redshift and thus to investigate the dark energy evolu- 
tion; P tp 2sz can be used to constrain the gas profile thanks to 
its dependence on the square of the gas density; Pk P sz, which 
we are going to discuss in more detail in this paper, is a tool 
to study the ICM dynamics and thus the history of the merging 
activity. 



Using a toy model, Diego et al. d2003l) showed that the 



kpSZ effect strongly depends on sub-clump motions, and that 
major merger events could make this signal equal or dominant 
on the other ones (depending on the assumed frequency). A 
different stud y based on a nume rical simulation has been car- 
ried out by lLavaux et alj d2004l) . who focussed their attention 
on the large scale structure of the kpSZ signal. Their results 
show that the environmental collapse into the potential of the 
forming cluster originates a characteristic kpSZ pattern with 
circular symmetry. They suggested to use this peculiar pattern 
to detect form ing clusters. In addition they confirmed the ana- 
lytic results of Sazonov & Sunvaev ( 1999b . but they did not in- 
vestigate the dependences of this signal on the merging activity 
of the clus ters and they d id not perform any extended statis- 
tical study. I Shimon et alJ (120061) . using a cluster simulated by 
the ENZO code, investigate the polarization levels and patterns 
produced by the different components originated by the kine- 
matic quadrupole moments induced by moving scattering elec- 
trons: their conclusion is that the signal is high enough (» lfiK) 
to become detectable in upcoming CMB experiments. 

Here we extend these previous works by analyzing two dif- 
ferent sets of high-resolution hydrodynamical N-body numeri- 
cal simulations, performed by assuming different modeling for 
the ICM physics, to investigate in detail the statistical proper- 
ties of the kpSZ signal. We demonstrate that it can vary up to 
one order of magnitude according to the internal cluster dy- 
namics. 

The plan of this paper is as follows. In Sect. [2] we intro- 
duce the Stokes parameters and the primary polarization of the 
CMB, then we present the main polarization contributions in- 
duced by galaxy clusters, discussing in more detail the kpSZ 
polarization. Sect.[3]is devoted to an extended description of the 
two sets of numerical hydrodynamical simulations used for our 
analysis. Sect. [4] introduces two different examples of galaxy 
clusters, a relaxed cluster participating to a major merger, and 
a dynamically more active object: the corresponding signals are 
compared to investigate the influence of the in-falling haloes on 
the kpSZ signal. The statistical analysis of the full samples of 
simulated objects is presented in Sect. [5] where we consider 
separately the results of non-radiative simulations and the ones 
obtained including different ICM modeling. Finally, our main 
conclusions are summarized in Sect. [6] 



2. Linear polarization by scattering with free 
electrons 

Scattering by free electrons induces a linear polarization 
component only if the incoming radiation has a quadrupole 
anisotropy. In the case of the CMB photons, we can distinguish 
different polarization components according to the origin of the 
quadrupole moment of the incoming radiation and to the loca- 
tion of the free electrons. In this section we briefly review these 
processes after having introduced the main quantities related to 
polarization. 



2.1. Parametrizing polarization 

Radiation can be parametrized by using the four Stokes pa- 
rameters: / represents the total intensity of the light, Q and U 
quantify the degree of linear polarization, while V gives the 
degree of circular polarization. Since the Thompson scattering 
does not introduce any circular polarization, we can set V = 
here. Notice that polarization is a pseudo-vector (i.e. it has not 
a verse), whose direction is given by 



1 

y = — tan 

2 



(1) 



The degree of polarization p can be expressed as 



Ve 2 + u 2 



(2) 



The Stokes parameters are additive along the line-of-sight, 
and thus the total effect can be computed as the integral of local 
effects: Q = / dQ and U = j dU. 

Throughout this paper, we compute the polarization effect 
assuming the incoming CMB radiation to be unaffected by the 
SZ effect, i.e. / =s Icmr'- this is justified by the assumption of 
low optical depth. 

2.2. Primary CMB polarization 

The CMB has a primordial polarization component originated 
at the last scattering surface. At that epoch, the universe was 
fully ionized and the quadrupole moment was related to the 
primary anisotropies. This signal is expected to peak inside 
the cosmological horizon scale at that time, because it is a 
causal process requiring a coherence into the electron motion. 
This polarization signal has been recently measured by the 
WMAP satellite, which also detected a large-scale component 
1(1 + 1)C™ 2 _ 6) /2tt = 0.086 + 0.029/iK 2 , interpreted as a signa- 
ture of a global re ionization occurred at redshift z = 10.9+H 



dPage et all 2006). 



2.3. Secondary CMB polarization by galaxy clusters 

Free electrons are also present into the ICM of galaxy clusters. 
In this case the polarization components depend on the nature 
of the quadrupole moment of the incoming radiation, as we are 
going to discuss in the following subsections. 
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2.3.1 . Polarization by single scattering of moving 
electrons: pkSZ 

Our work will be focused on the kpSZ effect which is discussed 
here in more detail. This effect arises from single scattering of 
the CMB photons by the ICM electrons. A moving electron 
sees the CMB radiation Doppler-shifted because of its peculiar 
motion and thus perceives a radiation intensity given by 



Cx 3 



,.vy(l+ft/) _ 1 ' 



(3) 



energy component. This polarization component depends on 
the optical depth t and on the amplitude of the quadrupole 
Q(z) at the cluster redshift, where the scattering process occurs: 
Pq « tQ(z). 

Because of that, if we have an estimate of t, for exam- 
ple thanks to SZ measurements, we can evaluate the intrin- 
sic quadrupole of the CMB at different redshifts. For a flat 
ACDM model with Oa = 0.7, the maximum signal is of the 
order of 4.9r //K. Assuming a typical value of t ~ 0.03, 
this leads to a maximum signal of ~ 0.14 /iK, as found 



where x = hv/(kT) is the adimensional frequency,/? = v e /c 
gives the electron velocity in units of the speed of light, y = 
(1 - y6 2 ) ^ 2 is the Lorentz factor, yu is the cosine of the angle 
between the direction of the electron velocity and the one of 
the incident CMB photon. The decomposition of this function 
into Legendre polynomials leads to 



/(v,Ai) = 



Cx 3 



1 



1\ e x (e x +l) , , 
I + hu + h \uf - - — -x 2 B 2 + 



■(4) 



The previous equation contains the quadrupole compo- 
nent necessary to produce a li near polarization orthogonal to 
the electron velocity direct ion (ISunvaev & Zerdovichi [1980a; 
Sazonov & SunvaevHl999l) . 

The resulting local contributions to the Q and U Stokes pa- 
rameters are 



dQ = -0.lo- T n e f(x)/3 2 cos(2x)dl 
dU = -Q.lo- T n e f(x)/3 2 sin(2x)dl 



(5) 



respectively. In the previous equation, cr T represents the 
Thompson cross-section, n e is the electron density, f3, is the 
component of jS on the sky plane (tangential velocity), x de- 
fines the position angle of the tangential velocity, and 



e*(e* + 1) 2 
fix) = r-X 

J 2(e x - l) 2 



(6) 



is the non-relativistic approximation of the frequency depen- 
dence of the polarization signal, as results of the quadrupole 
component of Eq.©. The Stokes parameters (see Eqj5j are nor- 
malized by the incoming mean intensity of the CMB radiation 
and are referred to the observer rest frame. They can be con- 
verted into a brightness temperature by dividing the previous 
quantities by 



*<*) = 



d ln/ v 



d In T e x - 1 



(7) 



analytically by ISazonov & Sunyaev (1999). This result has 
been confirmed by the numerical simulations performed by 
lAmblard & White! (120051) . A more detailed study of the primor- 
dial CMB quadrupole-induced polarization based on the com- 
bination of an analytic method w i th hy drodynamic simulations 
has been carried out by lLiu et al.l (120051) . who also investigated 
the contributions to the signal of different gas phases. Because 
of the dependence on Q, there are four orthogonal regions into 
the sky where Pq vanishes or at least is very low. These re- 
gions have been located by the WMAP satellite: one of them 
is at ga lactic coordinate (/, b) * (-80°, 60°), close to the Virgo 
cluster (ITegmark et all 12003). 

This observable is one of the most sensitive probe of the 
dark energy: in fact the growth of the gravitational potential, 
and thus of Q, is directly proportional to its equation of state. 
On the contrary all other related quantities depend at least on 
an integral of this quantity. 

2.3.3. Polarization by second scattering of the tSZ and 
kSZ quadruples: P p i tSZ , P p 2 ksz 

A first Thompson scattering of the CMB radiation with 
the ICM electrons introduces two secondary anisotropies, 
called kinetic and thermal Sunyaev Zel'dovich effects. They 
have a quadrupole component and thus a second scatter- 
ing would originate a linear polarization related to each of 
them. We call P p ^ t sz an d Pp 2 ksz tne polarization signals in- 
duced by the tSZ and kSZ effect s, respectively. As shown 
by ISunvaev & ZeFdovichl ( 1980a ). these polarization con- 
tributions depend only on the local density of free elec- 
trons and not on their temperature. The corresponding am- 
plitu des are proportional to th e square of the optical depth 
t 2 ( Sazonov & Sunyaev , 1999h . because of the two scattering 
events involved into the process: consequently they are very 
sensitive probes of the gas profile. The Pf-ksz is one order of 
magnitude stronger than the P p 2 tSZ , which can be even larger 
than the Pp signal, depending on the considered frequency 
dLavaux et all 120041) . 



Eq. (0 shows the kpSZ effect dependence on the square of the 
tangential velocity 1 and on its direction x- Therefore this ef- 
fect is very sensitive to the ICM dynamics and its amplitude 
can be used to investigate the merging activity. 

2.3.2. Polarization by single scattering of the intrinsic 
CMB quadrupole: P Q 

The CMB radiation has an intrinsic quadrupole moment given 
by the Sachs-Wolfe effect and by the late Integrated Sachs- 
Wolfe effect, which is mostly related to the presence of a dark separation through multi-band observations, similarly to the 



2.4. Possible contaminants 

Many other processes can introduce significant contributions 
to the linear polarization, like the Galactic synchrotron emis- 
sion by our Galaxy and by point radiosources, and the 
scattering with Galactic dust grains. In general these pro- 
cesses have different frequency dependences, allowing their 
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tech niques adopted for the CMB temperature maps (see, 
e.g., iTegmark et al. , 2000). The component-separation method 
could also be improved by including in the foreground sub- 
traction our knowledge on the distributions of dust, gas and 
magnetic fields in the Galaxy. In particular, if one is mostly in- 
terested in the polarization contributed by galaxy clusters, the 
observations could be focussed in regions at high galactic lati- 
tudes, where the foreground effect is minimal. Also the polar- 
ization originated by the CMB quadrupole {Pq) and the one 
induced by the thermal SZ effect {P p i t sz) have a different fre- 
quency dependence, favoring their separation. In addition Pq 
has typical scales larger than the kpSZ signal, so that the use of 
optimal spatial filters (see, e.g., lHuLl2001uMaturi et all 12005) 
could drastically reduce its contamination. 

The magnetic fields, combined with the presence of free 
electrons in galaxy clusters, may cause a Faraday rotation of 
the polarization orientation, reducing its importance. This ro- 
tation angle is relatively small also for massive g alaxy clus- 
ters: A RM = 1 - 10°(10G//z/v) 2 dOhno etalll2003l) . however it 
could become significant w hen considering our G alaxy at low 
galactic latitudes (see, e.g. jBurigana et aL , 20061) . 

It is well known that a polarization signal is also pro- 



duced by both the reionization process (see, e.g., Ng & Ng, 
1996 1 and weak gravita tional lensing effect (see, e.g., 
Zaldarriaga & Seliaklll998l) . In particular the deflection of the 



CMB photons leads to a smoothing of the small-scale features 
and to a power enhancement of the damping tail in the power 
spectrum. Notice that the analysis of the weak lensing effect on 
the CMB has be en suggested as possible tool to reconstruct the 
cluster masses jHuetall [2007): in this case the kpSZ signal 
represent a possible source of noise which must be taken under 
control to avoid systematic effects. 

In general all these contaminants should be analyzed in 
much more detail, but the promises of observing the kpSZ sig- 
nal in the next future are interesting. 

3. Numerical simulations 

For the purpose of this work we use two different sets of 
resimulated galaxy clusters. Both hav e been extracted fro m a 



dark matter(DM)-only simulation (see Yos hida et all 120011) of 
a "concordance" flat ACDM model, where the contributions to 
the density parameter from matter and cosmological constant 
are Q.o m = 0.3 and Qoa = 0.7 respectively; the Hubble constant 
in units of 100 km s -1 Mpc~' is h — 0.7. The initial conditions 
of the parent simulation, that followed the evolution of 512 3 
particles with mass of 6.86 x 10 10 hr x M & in a comoving cubic 
box of 479 hT x Mpc per side, were set considering a cold dark 
matter power spectrum normalized by assuming cr 8 = 0.9. 

In the catalogue of haloes identified in final output of the 
parent simulation by adopting a spherical overdensity criterion, 
we selected few objects to be re-simulated with higher mass 
and force resolution The new initial conditions have been cre- 
ated by using the "Z oomed Initial Conditions" (ZIC) technique 



(ITormen et all 11997). This method identifies in the initial do- 



main the Lagrangian region corresponding to each cluster up 
to few virial radii and populates it with more particles (of both 
DM and gas), properly adding small-scale modes. At the same 



time the volume outside the region of interest is sampled by 
low-resolution particles which allow to reproduce the original 
tidal effects of the cosmological environment. 

The new re-simulations, starting at redshift Zm ranging be- 
tween 35 and 60, have been carried out using for the first sam- 
ple the Tree-Smoot hed Particle Hydrodyn amics (SPH) paral- 
lel code GADGET dSpringel et all 1200 lh . while for the sec - 
ond set we ran its new version GADGET-2 dSpringeli 120051) . 
which includes an entropy- conserving formulation of SPH 
( Spri ngel & Hernquisn. [2002) and allow, by choice, the treat- 
ment of many physical processes affecting the baryonic com- 
ponent (see below). The gravitational softening length was kept 
fixed at 30 hr x kpc comoving (Plummer-equivalent) for z > 5, 
and was switched to a physical softening length of 5 hr x kpc at 
smaller redshifts. 

The first sample comprises 17 randomly selected galaxy 
clusters, whos e structure and d yna mical properties has been 



discussed in Ra sia et alJ (120041) and Tormen et al. (2004J), re 



spectively. These simulations, which follow only non-radiative 
hydrodynamics, assume a baryonic density corresponding to 
Q.[)b - 0.03; the particles masses in the high-resolution region 
range from 2 x 10 9 h~ l M Q to 6 x 10 9 h~ l M for DM, and from 
3 x 10 8 hT l Mv to 7 x 10 8 h~ l M Q for gas. The final objects have 
virial masses between 3.1 x 10 14 ^'M G and 1.7 x 10 15 h~ l M e , 
and virial r adifl between 1.4 and 2.5 h 1 Mpc. 

The second sample comprises 9 more isolated objects, 4 of 
them with the typical size of clusters (about 10 15 h~ l M Q ), and 
5 of group-like systems (about 10 14 /z _1 M o ). In this case the 
baryonic density has been fixed to Qo* = 0.04. Each cluster 
of this set has been re-simulated 4 times, turning on different 
kinds of physical processes: 

1 . gravitational heating only with an ordinary treatment of gas 
viscosity (hereafter ovisc model); 

2. like ovisc, but using an alternative implementation of the 
artificial viscosity (hereafter Ivisc model), which allows to 
better resolve the turbulence driven by fluid instabilities: 
this leads to a non-negligible contribution of non-thermal 
pressure support in the central region of galaxy clusters 



(Dolag et al., 2005); 



3. cooling, star formation and supernovae feedback with weak 
galactic winds, having a speed of ~ 340 km/s (csf model); 

4. like csf, but also i n cluding ther m al co nduction (see 



Jubel gas et ail 12004 iDolag et all |2004 . where the 



isotropic effective conductivity is fixed to 1/3 of the Spitzer 
rate (csfc model). 

A more detailed presentation of the properties of the clus- 
ters of this second set can be found elsewhere (see, e.g., 
Ettori et all l2006t iRoncarelli et alll2006l) . 

We want to stress here that the simulated clusters of the 
first sample have been randomly chosen from the complete list 
of haloes identified at z = in the parent simulation, with the 
only requirement of sampling in a sufficiently complete way 
the mass range above 10 14 /i _1 M Q . For this reason they have 
varying dynamical and environmental properties, comprising 



1 Virial radii are defined using the overdensity thre shold dictated by 



the spherical top-hat model (see, e.g., [Eke 



ensity thresho 
etall fl996). 
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relaxed and non-relaxed objects, isolated and interacting sys- 
tems. This allows us to consider this sample as representative 
of the whole population of galaxy clusters. On the other hand 
the clusters of the second set have been chosen to be more iso- 
lated structures, so that their main properties are not affected 
by recent major merger events and thus the effects of the four 
different ICM models can be better identified. 

To compute the final Stokes parameters related to the kpSZ 
signal, we need to integrate Eq. (O along the line-of-sight. This 
is done by considering boxes of (8 Mpc) 3 centered on the clus- 
ter: the resulting maps have a resolution of 15.6 kpc. The po- 
larization contributed by each particle, or moving cloud, adds 
constructively if their tangential velocities have the same di- 
rection, even if they have opposite versus. On the other hand, if 
their tangential velocities are orthogonal, the polarization tends 
to cancel. Because of that, only the coherent motion of gas 
clouds produces a polarization contribution, unlike the thermal 
motion of electrons, which is randomly oriented and does not 
contribute to the signal. 



4. The importance of the merging activity 

4.1. Two examples with different dynamical state 

In order to discuss the role of the merging activity, we ana- 
lyzed in more detail two clusters of our first sample, named 
g24+200 and g7. The first object represents a relaxed halo 
falling with high velocity into the potential well of another 
cluster. Consequently the velocity field of its gas is dominated 
by the bulk motion of the whole system. For this reason it can 
be considered at the same time like an example of a relaxed 
halo and of a major merger event because of its bulk motion 
given by the companion cluster. Conversely, the second object 
is not relaxed and is acquiring new matter from different in- 
falling subclumps having high velocities. 

Figure Q] shows a region of 8 x 8 Mpc 2 centered on clus- 
ter g7. This object has a virial mass of 1.46 x 10 15 hr x M & and 
a virial radius of 2.35 hr x Mpc (indicated by the circle in the 
plots). In the upper-left panel we overlay the mass-weighted 
tangential velocity field of the gas to its density iso-contours 
(red curves); the arrow length is proportional to the veloci- 
ties. The upper-right panel shows the corresponding map for 
the mass-weighted temperature of the gas; the color scale is 
given in the upper part of the panel. In the bottom left-panel we 
show the polarization pseudo-vector field, where the lengths of 
the lines are proportional to the polarization amplitude. Note 
that, for graphical purposes, the lines sample the signal on a 
grid having a lower resolution, giving only an idea of the polar- 
ization and velocity directions. Finally, we draw in the bottom- 
right panel the amplitude of the polarization signal using a log- 
arithmic scale (reported in the upper part). The polarization 
maps are computed at a frequency of 300 GHz. 

Since the polarization is proportional to the projected den- 
sity and to the squared tangential velocity (see Eq.|5]), the signal 
is peaked on the core of the merging haloes, where the col- 
umn density of electrons and the tangential bulk motion are 
large. The kpSZ signal of some in-falling blobs can be even 



stronger than the cluster core itself: an example is the blob ap- 
proximately located at the position (0.5 Mpc , 3.0 Mpc). 

In Figure [2] we display the same quantities for the cluster 
g24+200, which has a virial mass of 6.54 x lO u hr l M B and 
a virial radius of 1.80 h' 1 Mpc. Unlike cluster g7, for this ob- 
ject the polarization inside the virial radius is dominated by 
the cluster core where the signal shows a constant direction. 
Outside the virial radius there are two more evident structures 
which do not belong to this cluster: the first one is a small but 
fastly in-falling blob, the second one, located at the position 
(3.8 Mpc, 2.4 Mpc), represents a different cluster. The absolute 
value of the polarization signal of this cluster is large because 
it is taking part of a major merger event: in fact it is moving 
with a bulk velocity of 1200 km/s into the direction of the other 
cluster visible in this field. 

We discuss the differences between these two clusters and 
the importance of their merging activity into the following sub- 
sections. 

4.2. Excluding the merging activity 

In order to better understand the influence of the merging ac- 
tivity on the polarization signal, we compare the polarization 
map displayed in Figs.Q]and[2]with me ones obtained by artifi- 
cially "freezing" the internal dynamics. This is done by associ- 
ating the mean bulk velocity^ of the whole cluster to all parti- 
cles inside the box, so that the cluster moves like a rigid body. 
For simplicity hereafter we will refer to these simulations as 
"merging-free " simulations. 

The top panels in Figure[3]show the polarization signal ob- 
tained in such a way, for the unrelaxed cluster g7 (on the left) 
and for the relaxed cluster g24+200 (on the right). The po- 
larization signal is proportional to the length of the lines and 
it is displayed using the same scale adopted for Figs. Q] and 
[2] On the top of these panels, we also report an arrow corre- 
sponding to the bulk velocity of each galaxy cluster. The iso- 
polarization contours in the upper panels start from 1 nK and 
show levels in power of 2 steps, i.e 2, 4, 8, 16 nK, etc. Note that 
the large structure in the upper-right corner of the g24+200 
panel is an artifact due to the procedure performed to obtain 
the "merging-free" simulation. In fact, the bulk motion of the 
cluster could be assigned only to the particles inside the virial 
radius. Consequently, in the "merging-free" simulations only 
the region inside the virial radius has a physical meaning and 
will be considered in the following discussion. 

4.3. Quantifying the importance of the merging 

As described above, we computed the kpSZ signal of our 
two clusters, considering both for the real simulations and the 
"merging-free" versions of them. Comparing the correspond- 
ing results, we can evaluate the amplification of the polariza- 



2 In order to avoid the contamination given by in-falling blobs, this 
is defined as the velocity associated with the baricentre of the most 
inner 50 per cent of the cluster mass. 
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Fig. 1. Maps of a region of 8x8 Mpc 2 centered on the simulated cluster g7. Upper left panel: arrows represent the (mass-weighted) 
tangential velocity of the gas; the overlaid levels refer to the density iso-contours. Upper right panel: the gas (mass-weighted) 
temperature. Bottom left panel: the direction of the polarization pseudo-vector field. Bottom right panel: the polarized radiation 
amplitude. The different scale units or color levels are shown above the corresponding panels. In all panels, the circle indicates 
the virial radius. 



tion signal produced by the merging activity. This is done by 
computing the per cent increment of the signal, defined as 



A(0) = 



P(0) - P ml (0) 
P mi (0) 



(8) 



where P and P m f are the polarization signals given by the real 
and "merging-free" simulations, respectively. A map of this 
quantity is shown for both clusters in the bottom panels of 
Figure [3] Clearly the relaxed cluster in the left panel strongly 
differs from the non-relaxed object in the right panel. In the first 
case (g24+200) there is only a small increment of the signal (30 
per cent at most) and only over a small area centered on a small 
halo; the signal of the cluster core is substantially unchanged, 
showing an increment by only 10 per cent. In the second case 
(cluster g7) a large fraction of the object shows a substantial 



signal increment, especially in the outer regions, where most 
of the in-falling blobs are located. Here the merging activity 
produces an increment by 170 per cent in the cluster core and 
up to a factor of 10 at the positions of five different sub-halos, 
three of which have a signal even stronger than the cluster core 
itself. 

In order to better describe the distribution of the polariza- 
tion, in Figure|4]we report the fraction of the virial area covered 
by the galaxy clusters where the signal is larger than a given 
threshold. The results for the merging cluster g7 and for the 
relaxed system g24+200 are shown in the upper and bottom 
panels, respectively; in this case all three orthogonal projec- 
tions have been considered. The solid lines refer to the original 
simulations (i.e. including the merging activity), while the dot- 
ted lines correspond to the "merging-free" simulations. A more 
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Fig. 2. As Figure (03, but for the simulated galaxy cluster g24+200. 



quantitative evaluation of the differences between relaxed and 
un-relaxed clusters can be obtained by looking to the right pan- 
els, where we plot the fraction of the virial area as a function of 
the variation of the polarization signal produced by the merg- 
ing activity. It is evident that for the relaxed object g24+200 
there is not any significant polarization amplification from the 
inner dynamics: only 1 per cent of the cluster area shows an 
amplification of at least 30 per cent. On the contrary, the un- 
relaxed cluster g7 shows an increment of at least 300 per cent 
(800 per cent) on 10 per cent (1 per cent) of the cluster area; 
50 per cent of the cluster area inside the cluster has an ampli- 
fication of 50 per cent. Note that the distributions in the right 
panels do not reach the unity, because we considered only posi- 
tive variations. Thus the difference between unity and the max- 
imum value represents the fraction of the area within the virial 
radius with a negative A: even if it may be large, we checked 
that it corresponds to regions where the polarization signal is 
negligible. 

We notice that one of the three projections of cluster g7 
has a significantly lower bulk velocity with respect to the other 
ones (379 against 577 and 624km/s). As a consequence, in the 



Table 1. Per cent variation of the polarization signal for clusters 
g7 and g24+200 placed at redshift z — 0. 1 . The table reports the 
values corresponding to two different area fractions (10 and 1 
per cent) and to four different beam resolutions (ranging from 
to 10 arcmin). 
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merging-free simulation for this particular projection the po- 
larization signal is lower with respect to the other ones, but as 
high as the other cases when the merging activity is fully in- 
cluded. This shows that the gas dynamics originated during the 
merging events dominates over the cluster bulk motion. 

It is important to study how the previous results change 
considering different instrumental resolutions. This is evalu- 
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Fig. 3. Comparison with "merging-free" simulations. Left and right panels refer to the dynamically active cluster g7 and to the 
more relaxed cluster g24+200, respectively. Top panels: the polarization signal obtained from the "merging-free" simulations, 
shown by adopting the same scales and units of Figs. Q] and [2] The iso-polarization contours start from 1 nK and show levels 
power of 2 steps. An arrow proportional to the bulk velocity of each cluster and having its direction is displayed on the top right 
part. Bottom panels: maps for the per cent increment of the polarization signal, given by the merging activity, as expressed by 
Equation ((H). The corresponding color scale is shown above the panels. 



ated by placing the galaxy clusters at redshift z = 0.1 and con- 
volving the original maps with four different aperture beams. 
Table Q] summarizes the results obtained considering all three 
projections. Since the polarization signal contributed by the in- 
falling blobs is peaked on small scales, the convolution strongly 
smooths down the amplification produced by the merging ac- 
tivity. Obviously this effect is larger for lower resolutions, but 
it affects only the high-polarization region, leaving unchanged 
the signal on larger scales. 

Concluding this section, we notice that the example of these 
two clusters clearly shows the importance of the merging ac- 
tivity for the polarization signal and suggests two main ways 
to exploit this amplification to extract useful information. First, 
polarization could be used to investigate the ICM dynamics. 
In fact it could be possible in principle to obtain a complete 
three-dimensional reconstruction of the velocity field thanks to 



the combination with kSZ measurements (Na gai et al. . 2003 ): 
this would provide a picture of the merging history of our uni- 
verse. Second, the kpSZ signal could be used to determine 
the tangential bulk motion of clusters, allowing a determi- 
nation of the three-dimensional velocity power spectrum (al- 
ways in combination with kSZ data). This would lead to inde- 
pendent and co mplementary e s timates of cosmological param- 
eters (see, e.g.. [Tormen etall Il993l fMoscard ini et all [l996; 
Bhattach arva & Kosowskv , 2007b . In this case the merging ac- 
tivity represents a strong source of noise and a particular care 
should be used in selecting clusters. To discuss in more detail 
the reliability of these two possibilities, we need to repeat the 
previous analysis on more robust statistical bases. For this rea- 
son we analyze the entire samples of simulated clusters in the 
next section. 
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Fig. 4. The statistical distribution of the polarization signal for the merging cluster g7 (top panels) and for the relaxed cluster 
g24+200 (bottom panels). The results for three different cartesian projections of each simulated cluster are shown. The left panels 
present the fraction of the virial area covered by the cluster having a polarization signal larger than a given threshold (expressed 
in terms of AT/T); solid and dotted lines refer to results for the original simulations and for the 'merging-free' simulations, 
respectively. The right panels show the area fraction where the per cent variation of the polarization amplification produced by 
the merging activity (defined as in Equation[8]l is larger than a given threshold. 



5. Analysis of the full cluster samples 

5.1. Non-radiative simulations 

First, we perform a statistical analysis considering our first 
sample, composed by 17 different haloes simulated under the 
assumption of non-radiative hydrodynamics. For all objects, 
we consider the projections along three orthogonal lines-of- 
sights: this leads to 51 polarization maps. For simplicity we 
placed all systems at redshift z — 0.1. 

To investigate the effect of the instrument resolution, we 
convolve all polarization maps with a Gaussian kernel assum- 
ing different aperture beams. This is of primary importance be- 
cause the largest contributions to the total signal come from 
small haloes, which could be smeared out by a too large beam. 

In the left panel of Figure [5] we plot the cumulative dis- 
tribution of the kpSZ signal as a function of the covered area 
inside the virial radius. The different curves represent the me- 
dian of the 5 1 realizations and present the results obtained from 
the original maps (red line) and for different aperture beams: 
1' (green line), 5' (blue line) and 10' (black line). The thick 
lines refer to the simulations including the merging activity, 
while the thin lines are the results for the "merging-free" simu- 
lations, i.e. where only the cluster bulk motion is accounted for. 
Using the same line coding, in the right panel of Figure [5] we 
show the distribution of per cent amplifications for the polar- 
ization signal originated by the merging activity. As expected, 



the instrumental resolution causes a significant reduction of the 
maximum amplitude of the signal, decreasing also the incre- 
ment induced by merging. The minimum aperture beam affect- 
ing the signal is obviously related to the typical scale of the 
objects producing the polarization peaks. This scale is ~ 1 ar- 
cmin which, not surprisingly, is of the same order of the halo 
core radius (r c ~ 0.\\h~ l Mpc atz = 0.1). 

Our results suggest that the upcoming millimetric observa- 
tories, having angular resolutions of the order of FWHM 1' 
or better, will be sensitive to the merging activity of galaxy 
clusters, which is expected to amplify the kpSZ signal by a fac- 
tor of about 5 on at least 1 per cent of the virial area. The plot 
also shows that on average the cluster population has a max- 
imum signal of AT/T « 10~ 8 , while merging systems have 
peaks up to 10~ 7 : this confirms what we obtained from the anal- 
ysis of two clusters g7 and g24+200 in the previous section. 

We notice that this median distribution, which represents 
the typical expectation for the composite cluster population in- 
cluding both relaxed an d non-relaxed objec ts, is consistent with 
the results obtained by lLavaux et al.l ( 2004). However we point 
out that the merging clusters may have a kpSZ signal 10 times 
stronger. 

Additionally to the dynamic information carried on by this 
effect, it is important to notice that the kpSZ effect can be- 
come a non-negligible source of noise in the estimate of the 
CMB quadrupole when observing the polarization induced by 
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Fig. 5. The statistical distribution of the polarization signal as derived from the analysis of all projections of the 17 simulated 
galaxy clusters belonging to our first sample. Left panel: the fraction of the area having a polarization signal larger than a given 
threshold (expressed in terms of AT/T); thick and thin lines refer to results for the original simulations and for the 'merging-free' 
simulations, respectively. Right panel: the area fraction where the per cent amplification of the polarization signal produced by 
the merging activity is larger than a given threshold. Curves of different colors present the results for the original maps (red lines) 
and for different aperture beams: 1' (green line), 5' (blue line) and 10' (black line). 



the intrinsic CMB quadrupole Pq. This makes necessary both 
to have access to a component-separation technique based on 
multi-frequency observations and to carefully choose the clus- 
ters to be considered, which should be as relaxed as possible. 
This selection can be performed thanks to X-ray and SZ obser- 
vations. 

5.2. Effects of the ICM modeling 

Finally, we analyze the second simulation set to investigate the 
effects of different ICM modeling onto the statistical properties 
of the kpSZ signal. Again, the 9 objects have been projected 
onto 3 different orthogonal directions, leading to 27 polariza- 
tion maps for each of the ICM modeling described in Sect. [3] 

As an example, we shown in Figure|6]the polarization maps 
for the unrelaxed cluster g5 1 as obtained using the same pro- 
jection, but with the inclusion of different physical processes. 
This object has a mass of w 1 .3 x 10 15 h~ l M Q and a virial radius 
of » 2.3 hT l Mpc (corresponding to the circle in the plot). The 
properties of the large-scale polarization signal, coming from 
diffuse structures, are similar in the different panels, reaching 
values of the order of ~ Ink. We find significant changes only 
in the most central regions of the cluster, where the different 
ICM modeling plays an important role. In particular, from the 
comparison of the two upper panels, it is evident that the ovisc 
and Ivisc models give almost equivalent kpSZ signals, with 
consistent maximum amplitudes (967 nk and 1123 nk, respec- 
tively). Also the different features present in the map have co- 
incident positions and very similar appearence (they are only 
slightly smoother in the ovisc model). Since the only differ- 
ence in the ICM modeling implemented in these two simula- 
tions is the possibility of better resolving the turbulent motions 
driven by fluid instabilities in the Ivisc model, our results seem 
to suggest that their relevance for polarization is negligible. The 
processes of cooling, star formation and supernovae feedback 
with weak winds which are active in the csf simulation, pro- 
duce the effect of smearing out the highest peaks in the cluster 
core, strongly reducing the maximum signal (285 nk), almost 



half order of magnitude smaller than in the previous cases. The 
general properties of the polarization signal in the csfc model 
are similar to the csf case, with the exception at the small scales, 
where the signal is slightly smoothed down by the effect of the 
thermal conduction. 

In Figure|7]we show the polarization maps, but for the more 
relaxed and smaller cluster g6212 (having virial mass and ra- 
dius of « 1.15 x 10 14 /t'Mq and ~ 1.0 Mpc, respectively). 
In this case the kpSZ signal comes only from the most inter- 
nal regions (one third of virial radius) with peaks up to few 
tens of fjK. Again the maps for the two non-radiative simula- 
tions are very similar, with no significant difference between 
the ovisc and Ivisc models. For the simulations including cool- 
ing, star formation and feedback the typical values of AT/T at 
the centre are smaller. The reduction is particularly evident for 
the csfc model: unlike for the cluster g51, in the case of g6212, 
given its smaller mass, the thermal conduction is more efficient 
in smoothing down the signal. Comparing the results shown 
in Figs. [6] and [7] we find that, as expected, the different ICM 
modeling has larger effects on smaller clusters. However, this 
tendency is partially counterbalanced by the fact that the kpSZ 
signal is mostly coming from the cluster core and low-mass 
clusters are more concentrated than high-mass systems. 

A more quantitative picture can be drawn considering the 
statistical properties of total sample. In Figure [8] we show the 
distributions of polarization amplitude for each of the four ICM 
modeling here considered. In particular the curves represent the 
median of the 27 maps and refer to the differential and cumu- 
lative distributions (left and right panels, respectively). As pre- 
viously pointed out, the differences between the ovisc and Ivisc 
models are negligible. The radiative processes acting in the csf 
model produce some changes mostly in the high-signal tail, 
where the polarization amplitude is reduced by at most a factor 
of 2. A further small signal reduction is given by the switching 
on of the thermal conduction in the cfsc model. However, these 
differences are involving only few per cent of the cluster total 
area. 
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Fig. 6. Maps for the polarization signal for the same cluster, g51, but simulated with different ICM modeling: ovisc (top left), 
Ivisc (top right), csf (bottom left), csfc (bottom right). The same projection is shown in all panels; the circle represents the virial 
radius. 



6. Conclusions 

We studied the possibility of measuring the dynamics of 
the ICM during galaxy cluster mergers using the kinematic 
Sunyaev-ZeFdovich polarization (kpSZ) induced on the CMB 



1999). 



dSunvaev & Zel'dovichl Il980at ISazonov & Sunvaev , 
This effect takes origin from the single scattering of moving 
electrons: its dependence on the square of the tangential elec- 
tron velocity makes it a suitable probe to investigate the kine- 
matic properties of forming cosmic structures. 

We used the outputs of 53 high-resolution hydrodynami- 
cal N-body simulations to create realistic maps for this observ- 
able, which have been statistically analysed, focusing the atten- 
tion on the polarization signal induced by in-falling blobs. We 
demonstrated that the amplitude of the cluster signal can vary 
up to one order of magnitude depending on the internal cluster 
dynamics. 

We found that the median kpSZ signal of the whole pop- 
ulation is about 10~ & K, which could be observed by the up- 
coming millimetric observatories. In addition we showed that 



major merger events in unrelaxed clusters can produce a po- 
larization amplification of one order of magnitude, with values 
up to 10~ 7 K. These high peaks in the kpSZ signal are however 
related to fast in-falling blobs, corresponding to small angular 
scales: consequently, an high instrumental resolution (i.e. small 
aperture beam) is fundamental to avoid to smear out the signal. 

Our results show that the merging activity is of fundamental 
importance to estimate the three-dimensional velocity power 
spectrum of non-linear structure through the observations of 
the kpSZ and SZ signals. In fact, the amplification given by 
internal motions could be dominant with respect to the sig- 
nal originated by the bulk velocity of clusters, resulting in a 
non-negligible source of noise for such measurements. This is 
also true for the estimation of the CMB quadrupole from the 
observation of the polarization induced by the intrinsic CMB 
quadrupole Pq. This makes necessary both a robust compo- 
nent separation based on multi-frequency observations and a 
careful selection of the clusters to be observed, which should 
be as relaxed as possible. 
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Fig. 7. As Figure [6] but for the more relaxed and smaller cluster g6212. 
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Fig. 8. The differential (left panel) and cumulative (right panel) distributions of the polarization signal, as measured using all 
projections of the 9 simulated galaxy clusters belonging to our second sample. The different lines refer to results for simulations 
with different ICM modeling, as indicated in the panels. 
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